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Abstract 

Recent advances in photon detectors have resuhed in high-density imaging arrays that offer many performance 
and cost advantages. In particular, the excellent transit time spread of certain devices show promise to provide 
tangible benefits in applications such as Positron Emission Tomography (PET). Meanwhile, high-density, high- 
performance readout techniques have not kept on pace for exploiting these developments. 

Photodetector readout for next generation high event rate particle identification and time-resolved PET requires 
a highly-integrated, low-power, and cost-effective readout technique. We propose fast waveform sampling as a 
method that meets these criteria and demonstrate that sub-lOps resolution can be obtained for an existing device. 
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1. Background Summary of TOF System Timing Resolution 



Recent developments in high-density, high pre- 
cision timing photodetectors are finding appli- 
cations in Cherenkov detection techniques for 
particle identification, as well as medical imag- 
ing applications. This is true for both traditional 
vacuum-based [1] (Micro-Channel Plate, Hybrid 
Photo-Diode) and solid state [2] (Geiger-mode 
APD) photodetectors, and have resulted in lower 
cost, higher quantum efficiency, and improved 
transit-time-spread (TTS) photon detection op- 
tions. Fig. 1 illustrates that without these improve- 
ments, the performance of large-scale Time-Of- 
Flight (TOF) systems for precision spectrometers 
has stagnated at a resolution of about lOOps. 
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Fig. 1. A historical comparison of large Time-Of-Flight 
system timing performance, where progress has stagnated. 
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At the same time, the abihty for electronics to 
resolve at the 25ps or better level has been avail- 
able for decades, as reviewed in [15] and summa- 
rized in Table 1. 

Table 1 

An abridged summary, taken from [15], of precision time 
recording techniques that have been available for decades. 



Technique 


Resolution [ps] 


Citation 


TDC counter 


200 (20) 


[16] ([17]) 


Vernier (Chronotron) 


~ 15 


[18] 


Overlap Coin. (TAG) 


> 5 


[19] 


Start-Stop (TAG) 


> 5 


[20] 


Microwave 


< 1 


[21] 



For instance, consider a very simple Time to 
Amplitude Converter (TAG) circuit, as illustrated 
schematically in Fig. 2. When combined with a 
modest resolution ADC measurement, this is just 
about the simplest possible Time-to-Digital Con- 
verter (TDC). In its most reduced form the circuit 
consists simply of a current source, a switch, and 
a capacitor, as indicated in the figure. 

Vmax 

Vdd 




Fig. 2. Simple Time-to- Amplitude Converter circuit. 

At the start of some event of interest, the switch 
is closed and current starts to flow into the capac- 
itor, which is initially at OV. Because 



AQ = C-AV 



(1) 



and for constant current I, which is one ob- 



tains 



AV AT. 



(2) 



After the time interval AT, the switch opens and 
the voltage stored on the capacitor is a linear func- 
tion of the ratio of the current / over the capaci- 
tor C. Note that a maximum period Tmax can be 
measured for a maximum voltage Vmax. Rewrit- 
ing this expression we determine the time interval 
AT simply by measuring AV and having previ- 
ously determined the current and capacitance val- 
ues (or determining their ratio from calibration). 



AT=j-AV 



(3) 



Putting in numbers corresponding to typical val- 
ues for TOF applications and standard electronics, 
as listed in Table 2, one sees that 



SV 



Vmax 



= 1 mV. 



(4) 



Table 2 

Typical values for TOF TAC circuits using common elec- 
tronics components. 



Item 


Value 


Vmax 


4.095 V 


Capacitor 


10 pF 


Ramp Current 


409.5 M 


Full scale 


100 ns (typ.) 


ADC Resolution 


12-bits (4096 counts) 


Noise level 


ImV 



This is conveniently at about the amplitude 
often observed for board-level and other system 
pick-up noise in large systems. Therefore, the least 
count of this TAC is approximately 



6t = 



100 ns 

912 



25 ps 



(5) 



and for an ideal binary interpolation of this least 
count, one could do even a factor of — ^ better. 
Non-linearites due to imperfections and tempera- 
ture dependence in the current source and switch 
can degrade performance, but have been addressed 
in a number of commercial modules. For high rate 
environments, variants of this same basic scheme 
involving multiple ramps, have been employed for 
fast TOF measurement in an online trigger [22] , as 
well as hit logging prior to trigger formation [23] . 



2 



This performance is so much better than that 
for the ~ lOOps hmit mentioned earher that even if 
the TDC contribution were made zero, the overall 
timing would only improve by 3ps: 



0- = v/(100 ps)2 - (25 ps)2 97 ps 



(6) 



Therefore, the culprit must reside in the method 
of conveying the leading edge timing of the pho- 
todetector pulse into a signal that these proven 
methods can use. 



2. Amplitude Dependent Effects 

This well-known issue of "time walk" is illus- 
trated in Fig. 3, where a photodetector output sig- 
nal shape is amplitude dependent. Time walk cor- 
rects are used to compensate for the time of a 
threshold crossing slews to later times for smaller 
pulses. 
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Fig. 3. Example of the so-called time walk effect and its 
correction for the Belle TOF case. 

For detector signals with large amplitude varia- 
tions, the effect can be quite sizeable and must be 
corrected. A simple and proven method is to inte- 
grate the charge of the pulse within a gate around 
the signal of interest. While the leading edge of the 
signal can be complex, and can be modelled by a 
power-law, the observable of interest is the depen- 
dence on the pulse-integrated charge. This is illus- 
trated in Fig. 4, where the dependence observed 
for the Belle TOF system [11] is consistent with 



many previous observations as being proportional 
to approximately 
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Fig. 4. Amplitude dependence and result after applying 
TWC. 

Plotted as a function of ADC^^ in the upper 
left graph of Fig. 4, this dependence appears then 
as a linear expression. The effect of applying the 
correction is shown in the lower plots, where the 
amplude dependence has been removed and more 
than a factor of 4 in timing resolution improvement 
is observed. However, to further reduce scatter in 
the distribution requires a different or better pa- 
rameterization. 

For signals with smaller amplitude variations, 
such as in single photon detection. Constant Frac- 
tion Discrimination (CFD) techniques can be used 
instead, and may be combined with multi-hit de- 
vices [17,24] to achieve good timing performance. 
However, pile-up and rate effects need to be studied 
for each application. Moreover, all of the difficult 
work in making this circuitry functional is pushed 
onto the discriminator circuit. To reduce the effect 
of overdrive, the dependence of comparator out- 
put time upon amplitude above threshold, often a 
large amount of power may be required, making 
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the circuitry difficult to densely integrate. Due to 
the intrinsic memory of the CFD technique, it is 
susceptible to pile-up effects. 

A number of techniques, such as multi-level 
thresholding, have been proposed to directly mea- 
sure the risetime and provide a compensating cor- 
rection. Each threshold measurement may require 
an additional amplitude dependent correction to 
obtain the best possible timing. Adding more 
and more thresholding levels, one asymptotically 
approaches the penultimate technique, which is 
waveform sampling. 



Table 3 

Specifications of the BLABl ASIC used as a waveform sam- 
pling device to demonstrate the precision timing extraction 
technique. 



Item 


Value 


Sampling Input Channels 


1 


Storage rows 


128 


Storage cells/row 


512 


Total storage cells 


65,536 


Sampling speed (GSa/s) 


0.1 - 6.0 


Storage record 


10.9 - 655 Ats 


Operation mode 


continuous storage / readout 



3. Waveform Sampling 

A measurement complication at high event rates 
is pile-up, where the timing signal of interest may 
be superimposed on the tail of a preceding pulse. 
Without correction, and often being unaware that 
this superposition is present in the data, an incor- 
rect timing determination will be made. If afford- 
able and fast enough [25] , and if the aperture jitter 
is small enough, waveform sampling provides max- 
imum information extraction. The benefits may be 
summarized as follows: 

- Pile-up impact extraction 

- Fit to actual pulse shape 

- Average to reduce single measurement errors 

- Observe changes to pulse shape over time 

where the last issue can be particularly impor- 
tant for detectors that undergo radiation damage 
or aging effects such that the intrinsic pulse shape 
evolves with time. 



4. The BLABl Monolithic Sampler 

To demonstrate the viability of the waveform 
sampling technique, we use the first Buffered 
LABRADOR (BLABl) ASIC [25]. This is a 2^6 
(64k) cell deep waveform sampling ASIC that was 
developed as a deep-sampling extension of the 
successful LABRADOR ASIC [26]. The key pa- 
rameters of this device are summarized in Table 3. 



The evaluation system consists of the small cir- 
cuitboard, containing 2 BLABl ASICs, and read 
out over USB2.0 by a laptop running a custom 
data acquisition program based on the wx Widgets 
tool kit. A photo of this compact measurement sys- 
tem is displayed in Fig. 5. This compact configu- 
ration can turn any PC (or laptop) into a high- 
performance digital signal oscilloscope. 




r 



Fig. 5. Photograph of the "oscilloscope on a chip" test 
set-up, where the 2 BLABl ASICs are at the left, near the 
SMA input connectors. 
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5. Precision Timing Performance 

All of the standard precision timing measure- 
ment techniques discussed at the outset suffer from 
practical limitations in actual application, which 
has served to degrade the realized large system per- 
formance. In the end, one simply cannot do better 
than having a high-fidelity "oscilloscope on a chip" 
for every sensor channel. Cost and data volume 
precluded this type of waveform recording until 
recent generations of SCA ASICs [26,27,28,29,30] 
demonstrated such techniques were practical, es- 
pecially for large systems. 

5.1. Calibration 

In order to achieve fine timing resolution, a num- 
ber of calibrations need to be performed and cor- 
rections need to be applied. As reported in ear- 
lier devices, the temperature dependence [25,26] of 
the sampling baseline must be compensated. The 
effect of this 0.2% /°C temperature dependence is 
seen in Fig. 6. 
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Fig. 6. Event-by-event sampling speed comparisons of two 
BLABl ASICs indicating a common temperature depen- 
dence term. The sampling speed is continually monitored 
and used in calibration to achieve the maximum timing 
performance. 

In order to address bin-by-bin timing width 
differences, a couple of different calibration tech- 
niques have been tried. The first utilizes a sine 
wave zero-crossing technique used for calibrating 
the LABS ASIC [26]. That technique works best 



when the frequency of the sine wave is such that 
the measured interval between zero crossings can 
be uniquely assigned to a limited number of bins 
between successive crossings. 




Fig. 7. Aperture width determination for individual sample 
cells via the histogram occupancy technique described in 
the text. 

Due to sine wave curvature, this technique has 
an irreducible systematic error that is a function 
of sample rate. A more successful technique is to 
histogram the zero crossings of a sine wave and use 
the bin occupancy to derive the effective aperture 
width, as may be seen in Fig. 7. 

When a sampling strobe is inserted into the 
BLABl, the strobe propagates across a Switched 
Capacitor Array (SCA) row and exits the ASIC via 
a monitor pin. The sampling speed is determined 
by measured this propagation delay, a method 
implemented by creating a feedback loop with the 
insertion-extraction chain and an inverter located 
inside the companion FPGA. This feedback in- 
verter is only connected after waveform readout of 
a triggered event and then disconnected once the 
sampling speed has been measured for that par- 
ticular triggered event. A so-called Ripple Carry 
Out (RCO) period may then be expressed as 



^total 2 * (^latency -^waveform) 



(7) 



where the factor of 2 comes from the FGPA feed- 
back inverter toggling the RCO period, T„avoform 
is the SCA propagation time, and Tiatcncy is the 
additional time accrued due to internal FPGA 
routing, circuit board routing, and routing to and 
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from the SCA within the ASIC. To measure Ttotai 
an FPGA-based Time-to-Digital Converter (TDC) 
is used [31]. Both phases of a 250 MHz clock are 
used by the TDC to obtain a least count of 2ns. 
Improved resolution is obtained by integrating 
the RCO period for 2048 cycles, corresponding to 
an aggragate single-cycle least count of 0.977 ps. 
Tlatoncy IS obtained by subtracting Twavcform from 
Tlotai, where ^waveform IS determined from a fit to 
a fixed frequency sinusoidal signal. The sampling 
speed is then 



/SPD = (T'v 



waveform 



/RoWiongth) ^ 



(8) 



where /spd is the sampling speed and Rowiength 
is the number of SCA cells in a row, which for the 
BLABl ASIC is 512 cells. As reported in Ref. [25], 
the sampling timebase is temperature dependent. 
This effect may be seen as the correlated portion of 
Fig. 6. By measuring the BLABl timebase, these 
timebase shifts can be monitored and corrected of- 
fline. 

5.2. Ideal Test Signals 

An idealized test pulse of fixed amplitude and 
width was generated by taking an Agilent 33250A 
arbitrary waveform generator pulse output, AC 
coupling the signal with an inline capacitor, RF 
splitting the signal, and waveform sampling both 
signals using two separate BLABl ASICs with ad- 
ditional cable delay for the second channel. This 
AC coupling served as a filter and shaped the pulse 
into an almost gaussian shape, as seen in Fig. 8. 

The result obtained for taking the time differ- 
ence of a simple gaussian fit to the signals in both 
channels is plotted in Fig. 9. 

5.3. Varying Sampling Window Position 

The two pulses were then recorded without re- 
gard to the sampling strobe phase. This allows 
the two pulses to move around within the wave- 
form sampling window. During the analysis, a bin- 
by-bin timebase correction and an event-by-event 
sampling speed correction were applied. The time 
was measured by taking the cross correlation of the 
two waveforms. The cross correlation was done by 
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Fig. 8. Pulses recorded in two separate BLABl devices, 
with a small timing offset. A gaussian fit, as shown in the 
bottom waveform, is used to extract the event-by-event 
time of each signal. 



taking the Fast Fourier Transform (FFT) of both 
waveforms, complex multiplying one FFT by the 
complex conjugate of the other FTT, and Inverse 
Fast Fourier Transforming (IFFT) this product to 
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time (ns) 

Fig. 9. Time difference between a pair of idealized pulses. 
Each edge can be inferred to bo extracted a factor of \/2 
better, or about 4.5ps RMS. 

get the cross correlation waveform in the time do- 
main. The time position of the maximum point 
in the cross correlation waveform is the time dif- 
ference between the best estimator of the location 
of the two pulses. Because the sampling speed be- 
tween the two ASIC chips vary from each other and 
the time spacing between the bins is not constant, 
before FFT a fifth order data spline is applied at 
a fixed time spacing very close to the average time 
interval. This is done to keep the sampling spacing 
between points constant and minimize frequency 
domain distortion. To get finer time spacing in 
the cross correlation waveform without distortion, 
we symmetrically insert zeros in the cross correla- 
tion FFT for frequencies greater than the Nyquist 
frequency and less than the negative Nyquist fre- 
quency before the IFFT operation. The time res- 
olution obtained using this technique is shown in 
Fig. 10, corresponding to an RMS of 6.7 RMS or a 
a of 6.0 ps. 

For complex curvature along the leading edge of 
the signal, the timing resolution obtained is seen to 
be rather sensitive to the method choosen to char- 
acterize the signal "hit" time. Unless the photode- 
tector signal is for a single p.e. quanta, the actual 
shape can be rather complex and dependent upon 
photon arrival statistics. Even in this simple case, 
noise and aperture systematics upon the leading 
edge can be important and can also be reduced 




time (ns) 



Fig. 10. Measured time resolution by taking the cross cor- 
relation of two ideal pulse waveforms. 

by using multiple samples to fit to an analytic sig- 
nal shape. In general, the estimate error can im- 
prove as something like 1/ViV for N samples along 
the leading edge. This is perhaps the most pow- 
erful aspect of having the full waveform samples 
to fit. Individual sampling errors can be averaged 
out. Examples are provided in the following sub- 
section, where it is clear that at the sampling rates 
being studied, this waveform recording technique 
logs many samples on the leading edge, which can 
be used to improve the signal timing extraction. 

5.4. PMT signal observation 

A convenient feature of the BLABl ASIC is that 
a PMT output transmitted over a 50ri coaxial ca- 
ble can be directly connected to the BLABl input. 
A waveform recorded from a Burle 85011 Micro- 
Channel Plate PMT is displayed in Fig. 11. 

This near ideal waveform is utilized as the simu- 
lation template in the following section on perfor- 
mance limits from simulation. 

5.5. MCP-PMT signal timing 

A set of 10k triggered (cosmic and dark count) 
events were recorded, as illustrated for an exam- 
ple signal in Fig. 11. The signal amplitude was al- 
lowed to fluctuate and the measured distribution 
of integrated charge is shown in Fig. 12. 
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40 42.5 45 
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Fig. 11. Example waveform used in this analysis as recorded 
from a Burle 85011 Micro-Channel Plate PMT. 
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Fig. 12. Pulse height distribution for the signals used in 
this timing analysis. 

Ideally one would use a fast timing (ps) laser to 
do this measurement, however the Transit Time 
Spread of the photodetector would then be an is- 
sue. Instead we use the same pulse, split it, and 
record the signal on two storage channels. After ap- 
plying timebase and time bin corrections, the raw 
timing difference distribution seen in Fig. 13 is ob- 
tained. 

This distribution is rather wide and is related to 
a residual timing miscalibration. This can be seen 
by looking at the time difference as a function of 
the peak time determined on the first channel, as 
shown in Fig. 14. 

We created a secondary bin-by-bin timing cor- 
rection related to the time position of the first peak 
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Fig. 13. Raw timing performance after simple timebase 
correction. 
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Fig. 14. Observed residual timing structure, indicating a 
limitation in the previous calibration procedure. 

using the data from Fig. 14. After applying this sec- 
ondary tinming correction, Fig. 15 illustrates the 
improvement in timing performance. 

For this measurement, the amplitude of the 
MCP signal varies from event to event. For the 
waveform data, we calculate the charge of the 
signal by summing all the waveform data points, 
multiplying by the average timebase, and divid- 
ing by the termination resisitor. After applying 
all the corrections previously mentioned, a profile 
histogram of the charge versus time in Fig. 15 
is shown in Fig. 16, which illustrates that only 
a slight amount of residual time walk is remains 
after applying the correction. 
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Fig. 15. Time distribution after applying an optimized, 
binned timing correction. 
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Fig. 16. Improved time difference residual as a function of 
integrated pulse charge. 

6. Limitations 



While the analog bandwidth of the BLABl is 
adequate for many RF recording applications, a 
higher bandwidth device will be explored, based 
upon the lessons learned from this first device. In 
particular, the fanout structure and design of the 
analog amplifier tree is being scrutinized and im- 
proved in simulation. It is hoped that an almost 
arbitrarily large storage depth can be accommo- 
dated up to IGHz of analog bandwidth through a 
careful layout of the buffer amplifier cascade array. 

To gain insight into the parameters that are most 
important to attaining the highest possible tim- 



ing precision, a series of simulation studies were 
performed. In each run a series of 10k events were 
generated with a time offset equal to that used in 
the previous section. A real pulse template shown 
in Fig. 11 has been used, and is perturbed as de- 
scribed for each of the following cases. The simula- 
tion outputs the distorted waveforms into the same 
data structure as used by the acquisition software 
and the same cross-correlation method described 
earlier is applied. 

6.1. Random Noise 

The result due to the random addition of noise 
to each sample of the two waveforms is illustrated 
in Fig. 17. Comparing the curve with the measured 
result, an effective noise of ~ 3mV is indicated. 

Siinulated Time Difference Dependence on Noise 




2 3 4 

Added RMS noise [mV] 

Fig. 17. Fitted time resolution as a function of added noise. 

However, since we know that the noise in each 
sample is only about ImV, additional contribu- 
tions need to be considered. 



6.2. Sample Aperture Jitter 

Another expected contribution to the extraction 
of the timing resolution is jitter on the timing sam- 
ples. Using the same simulation configuration, the 
sample noise is now fixed and the effect of adding 
a randomized aperture jitter is plotted in Fig. 18. 

The local signal slope is used to determine a 
voltage offset expected by sampling at the wrong 
time. Since the sampling period is small compared 
with the typical voltage step between samples, a 
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Simulated Time Difference Dependence 
(1mV RiVIS noise) 



6.3. Sampling Timebase Jitter 
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Fig. 18. Timing resolution with ImV noise as a function 
of aperture jitter. 



rather large jitter is needed to reproduce the ob- 
served time resolution. From direct measurements 
of timebase jitter, we expect a far smaller value 
than is required to explain the time resolution ob- 
served. It is noted that the two effects are degener- 
ate in the sense that they lead to a sample voltage 
dispersion from the ideal value. Considering this, 
one can express the time difference resolution in 
terms of a combined value characterized in terms 
of a Signal to Noise Ratio (SNR). This relationship 
is seen in Fig. 19, where it is clearly seen that an 
effective SNR of 100 or greater is desired to obtain 
time resolutions of 5ps or better. 



Time Difference Dependence on Signai-Noise Ratio (SNR) 



100 
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Fig. 19. Impact of effective sampling noise on timing reso- 
lution, expressed in SNR. 



While the pedagogical illustrations above are 
helpful in understanding what limits the time res- 
olution that can be obtained, there is an additional 
contribution that is an artifact of the use of the 
BLABl ASIC. Specifically, it contains no dedicated 
Delay Lock Loop (DLL) circuitry, and thus the 
sampling timebase is sensitive at the level of about 
0.2%/°C [25]. Even though an event-by-event cor- 
rection is performed, as described previously, it is 
quite probable that a residual correction error con- 
tributes to the time resolution. A simulation is per- 
formed as shown in Fig. 20, where ImV of noise is 
added in quadrature with Ips of sample aperture 
jitter on top of an overall event-by-event sampling 
timebase jitter. This jitter is assumed to be uncor- 
related between the two ASICs. 



Simulated Time Difference Dependence 
(timebase jitter study) 




0.005 0.01 0.015 0.02 0.025 

Baseline jitter [%] 

Fig. 20. Time resolution with ImV noise, Ips of sample jit- 
ter, and as a function of event-by-event uncorrected time- 
base jitter. 



In the limit that the known ImV noise only con- 
tributes 2ps to the overall time resolution and that 
the sample aperture jitter is negligible, this effect 
is seen to be dominant. Quantitatively, this cor- 
responds to about a 10% error on correcting the 
event-by-event sampling timebase. If true, this can 
be improved through better timebase stabilization 
in a next version of the ASIC [25], or through the 
use of other devices as described next. 
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7. Other Implementations 



The waveform sampling technique described is 
more general than the particular ASIC presented. 
Indeed if price and power arc no object, a num- 
ber of commercial digitizers are available, with new 
devices continually coming to market. In general 
those devices, besides being rather power himgry, 
also suffer from a limited dynamic range. However, 
for many applications, they are already or will be- 
come viable. 

In the near term, a number of custom ASICs 
with better analog resolution [32], faster sam- 
pling [33] and faster readout [34] , will be available 
to implement these waveform sampling tech- 
niques. Further studies of the performance limits 
on timing signal extraction using this waveform 
sampling technique is a high priority for the follow- 
on BLAB2 ASIC [25] as well. In the longer term, 
these low-cost devices show great promise to re- 
place standard ADC and TDC techniques in many 
future high energy and particle astrophysics detec- 
tors [35]. As we have demonstrated, they should 
be able to do so with excellent timing precision. 



8. Summary 

A waveform sampling technique has been shown 
to obtain sub-lOps timing resolution for readout of 
photodetector signals. Measurements were made 
with an existing compact, low-power and pipelined 
readout ASIC. Such a technique can be scaled up 
to very large channel counts in a cost-effective way. 
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